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Abstract: SARAH is a Mathematica package for building and studying supersym- 
metric models. It calculates for a given superpotential and gauge sector the full 
Lagrangian of a model. With the new version of SARAH it is possible to calculate 
automatically all interactions for the different eigenstates and write model files for 
FeynArts and CompHep/CalcHep. In addition, the tadpole equations are calculated, 
gauge fixing terms can be given and ghost interactions are added, particles can be 
integrated out and non supersymmetric limits of the theory can be chosen. CP and 
flavor violation can easily be switched on or off. 



1. Introduction 



Supersymmetry is one of the most popular extensions of the standard model (SM) 
of particle physics [p], §, §, §, [|, §, 0] . The minimal supersymmetric standard model 
(MSSM) is well studied and there are several tools to explore the model: For every 
event generator or diagram calculator exists a model file, thus the MSSM can be 
handled out of the box. But what happens, if you want to change the model, extend 
the particle content, add a new gauge group or drop some assumptions about the 
parameters? Before you are able to do really physics and phenomenology you have 
to do a lot of work before: Check if your new model is free from gauge anomalies, get 
the Lagrangian out of the superpotential and extract all interactions, calculate the 
mass matrices and diagonalize them to get your mass spectrum, derive the tadpole 
equations to find the minimum conditions for your potential and so on. After all, 
if you are content with the first impression of your model, you might want to do 
phenomenological studies. But before you can use one of the existing programs, you 
have first to write a model file for them, what is again a very time consuming work. 
This is exactly that kind of work SARAH can now do for you. SARAH just needs 
the gauge structure, particle content and superpotential to produce all information 
about the gauge eigenstates. If you have also broken gauge symmetries and particle 
mixings you can easily add them. Also the gauge fixing terms can be specified. At 
the end you have the choice, if you want a F/TjrjX file with all information about 
your model, if a model file for FeynArts |13| or CalcHep/CompHep [[L4| should be 
written or if you just want to check some things using Mathematica. This model 
files can of course also be used with FormCalc [|ll[] for calculating Feynman diagrams 
with Mathematica or micrOmegas [0[ for dark matter studies. The intention by 
the development of SARAH was to make it very flexible: There is a big freedom 
for the matter and gauge sector you can handle. The work with SARAH should be 
easy: Every information SARAH needs are specified in an easy to modify model file. 
Nevertheless, SARAH is also fast: A existing model can be changed within minutes, 
and the needed time for doing all necessary calculations and writing a model file is 
normally less than 10 minutes. 



2. Overview 

SARAH is a package for Mathematica 1 and was tested with versions 5.2 and 7.0. 
The package can be downloaded from 

http : / /theorie . physik . uni-wuerzburg . de/~f nstaub/ sarah . html 

1 Mathematica is a protected product by Wolfram Research 
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and a comprehensive manual for SARAH can be found here ||. 

After you have extracted and saved the package, you can load SARAH in Math- 
ematica by 

«"\Package Directory\SARAH.m" 

and start the calculations for a model with 
Start ["modelname"] ; 

The following steps are done automatically: 

• The model is checked for gauge anomalies 

• The Lagrangian for the gauge and mass eigenstates is calculated 

• Ghost terms are calculated from the gauge fixing terms 

• If chosen, particles are integrated out 

• Mass matrices are derived 

• The tadpole equations are calculated 

• All tree level masses are calculated 

• A spectrum file is read in and the mixing matrices are calculated 

On modern computers, this will last only a few minutes. After the automatized 
calculations are finished, you can just play with the model by checking distinct 
vertices, look at the different mass matrices, give numerical values to the parameters 
and calculate the mass spectrum, or you can generate model files for FeynArts and 
CompHep/CalcHep for further phenomenological studies. Also a FTgX output with 
the particle content, mass matrices, tadpole equations and all interactions including 
corresponding Feynman diagrams can be written. 

3. What Models are possible 

In principle, SARAH can handle every N = 1-SUSY theory, in which the chiral 
superfields are arranged in SU(2) doublets and singlets. The gauge sector can be 
extended by abelian and also non abelian gauge groups, more particles and/or gener- 
ations of particles can be added to the matter content, one or more gauge symmetry 
breakings can be defined, particles can be integrated out to get an effective theory or 
just 'deleted' to get a non supersymmetric limit of a theory like the standard model. 
Non canonical terms can be added to the Lagrangians and existing interactions can 
be redefined. 



2 



You have always full control over the properties of the parameters of your model: 
You can switch on and off flavor and CP violation independently for all parameters, 
define relation between the parameters, add numerical values or specify the position 
in a Les Houches input file. 

In the package archive of SARAH already some models are included: The MSSM, the 
MSSM in CKM basis, the NMSSM, the /xi/SSM, the MSSM with bilinear R-parity 
violation, the SM, an effective MSSM with integrated out Gluino and a MSSM with 
additional U(l). 

4. The Model Files 

All information about a model are saved in three different files: Model .m, parameters .m 
and particles .m. Only the first one is absolute necessary and contains the infor- 
mation about the gauge sector, particle content, superpotential and mixings. In 
parameters .m you can assign properties to all parameters of your model, give nu- 
merically values and define the PTj^X name for each parameter. In particles. m 
additional information about the particles are given, which might be needed for an 
appropriate output: R-parity, mass, width, PDG, PTjtX name and output name. 
All three files are written in an intuitive way and can easily and quickly be changed. 
This, we want to demonstrate at an example in the appendix. 

5. The Output 

With 

ModelOutput[EWSB,WriteTeX->True, WriteCHep->True , WriteFeynArts->True] ; 

it is possible to create model files for FeynArts and CalcHep/CompHep as well as a 
P/T^X file with all information about the model for given eigenstates. 

5.1 WVyX 

In the generated PTjtX files are all information about a model for one set of eigen- 
states included: Particle content, mixing matrices, tadpole equations and all in- 
teractions. For the interactions are also automatically the corresponding Feynman 
diagrams are drawn by using the PTgX package FeynMF of Thorsten Ohl flnj. For 
generating the pdf file from the different PTjtX files and all Feynman diagrams a 
batch file is written automatically by SARAH for Linux and for Windows. 

5.2 CalcHep/CompHep 

There have to be taken some difficulties into account when writing a model file for 
CompHep or CalcHep: The color structure is in implicit, therefore a four point in- 
teraction of four colored particles can't be handled. This was solved by the authors 
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of CalcHep/CompHep by adding auxiliary fields and splitting such interactions in 
two parts. SARAH does the same splitting by taking the results from the F- and 
D- Terms. This has also the nice advantage that not all possible combinations of 4 
scalar have to be written separately and leads to an increasing of the speed writing 
the model file. 

In addition, CalcHep/CompHep works with three different kind of ghosts: Fadeev 
Popov ghost, Goldstone ghosts and tensor ghosts. The first two are supported by 
SARAH and the vertices are calculated dynamically, i.e. the model files produced by 
SARAH support both Feynman gauge as well as unitarity gauge. The tensor ghost 
of the gluino and its one interaction is hard coded. 

CalcHep/CompHep doesn't support complex parameters. Therefore SARAH has 
an option to split each vertex in a real and imaginary part in the Lagrange file of 
CalcHep/ CompHep in order to enable also theories with CP violation. 



5.3 Fey n Arts 

Generating the FeynArts model file is more straight forward than the files for CalcHep/ 
CompHep. But in addition to the normal model file for FeynArts SARAH writes 
a second file including additional information about the model: All defined depen- 
dences, the numerical values for the parameters and masses if they are available and 
special abbreviations to speed up the calculations with FormCalc. 



6. Checks 

We have checked the model files generated with SARAH for the MSSM against the 
existing files of FeynArts and CalcHep both at the vertex level and at the process 
level. This means, we have compared the value of each vertex for a given set of 
parameters and all possible combinations of generations (more than 15000), and we 
have also calculated several processes with the old and new model files. In addition we 
have calculated the relic density with micrOmegas for two different set of parameters. 
A recent version of SARAH was also used for PI and the results were cross checked 



in several ways. Also one of the authors of |L6] checked the results of SARAH for 
the [iv SSM. 
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A. Example: From MSSM via NMSSM to /iz/SSM and more 

A.l The MSSM 
A.l.l The Model File 

First, we want to explain the different parts of the model file of the MSSM. 

1. The gauge sector is U(l) x SU{2) x SU(3) and is just defined by adding the 
corresponding vector superfields. 

Gauge [[1]]={B, 1, hypercharge, gl}; 
Gauge[[2]]={WB, 2, left, g2}; 
Gauge [ [3] ] ={G , 3, color, g3}; 

The different parts are: The superfield name, the dimension of the gauge group, 
the name of the gauge group and the name of the gauge coupling. 

2. The doublet superfields are q, I, Hj, H u and added by 

Doub[[l]] = {uL, dL, 3, q, 1/6, 1, 1}; 
Doub[[2]] = {vL, eL, 3, 1, -1/2, 1, 0}; 
Doub[[3]] = {HdO, Hdm, 1, Hd, -1/2, 1, 0}; 
Doub[[4]] = {Hup, HuO, 1, Hu, 1/2, 1, 0}; 

The different parts are: The name of the up component, the name of the down 
component, the number of generations, the name of the superfield the different 
charges under the gauge groups. 

3. The singlets superfields are d, u, e and added by 

Sing[[l]] = {dR, 3, d, 1/3, 0, -1}; 
Sing [[2]] = {uR, 3, u, -2/3, 0, -1}; 
Sing [[3]] = {eR, 3, e, 1, 0, 0}; 

The definition is analog to the definition of the doublets. 

4. The trilinear terms of the superpotential leading to the Yukawa interactions 
are 

W tri = -LY e eH d ~ QY d dH d + QY u uH u 
and given in SARAH by 
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TriW = {{q, Hu, u, Yu,l}, {q, Hd, d, Yd, -1}, {1, Hd, e, Ye, -1}}; 

5. The bilinear term leading to the Higgsino mass is 

W bi = fiH d H u 

This reads in SARAH 

BiW = {{Hu,Hd,\[Mu] ,1}}; 

6. The gauge fixing terms for the unbroken gauge groups are 

GaugeFix={ {Der[VWB] 5 -1/(2 RXi[VWB])}, 
{Der[VG] , -1/(2 RXi [VG] ) }}; 

This has the meaning of 

7. Begin with the definition for the first set of rotations 
Rotation [[1]]={ 

(a) The name of the rotation is 
EWSB 

(b) The vector bosons and gauginos rotate as follows 

{ {VWB, {l,{VWm,l/Sqrt[2]},{conj [VWm] ,1/Sqrt [2]}}, 
{2 , {VWm , -\ [Imaginaryl] /Sqrt [2] } , 

{conj [VWm] , \ [Imaginaryl] /Sqrt [2] }} , 
{3, {VP, Sin[ThetaW]},{VZ, Cos [ThetaW] }}} , 
{VB, {1,{VP, Cos [ThetaW] },{VZ, -Sin [ThetaW]}}}, 
{fWB, {l,{fWm,l/Sqrt[2]}, {f Wp, 1/Sqrt [2] }} , 
{2 , {f Wm , -\ [Imaginaryl] /Sqrt [2] } , 
{f Wp , \ [Imaginaryl] /Sqrt [2] }} , 
{3,{fW0,l}}} 

}, 

This is the common mixing of vector bosons and gauginos after EWSB: 



W ± = -^Wi =F iW 2 , W ± = -^Wi =F iW 2 
v2 v 2 

Z = — sin Q W B + cos <d w W 3 , A = sin <dwW 3 + cos <d w B 
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(c) The neutral components of the scalar Higgs get Vacuum Expectation Val- 
ues (VEVs) 

H d 

K _ 

Here v u and v d are the VEVs of the Higgs. This is added to SARAH by 

{{SHdO, {vd, 1/Sqrt[2]}, {sigmad, \ [Imaginaryl] /Sqrt [2] } , 

{phid, 1/Sqrt[2]}}, 
{SHuO, {vu, 1/Sqrt[2]}, {sigmau, \ [Imaginaryl] /Sqrt [2] }, 

{phiu, 1/Sqrt[2]}}} 

(d) The particles mix after EWSB to new mass eigenstates 

{{{SdL, SdR}, {Sd, ZD}}, 
{{SuL, SuR}, {Su, ZU}}, 
{{SeL, SeR}, {Se, ZE}}, 
{{phiu, phid}, {hh, ZH}}, 
{{sigmau, sigmad}, {Ah, ZA}}, 
{{conj [SHup] ,SHdm},{Hpm,ZP}}, 
{{fB, fWO, FHdO, FHuO}, {LO, ZN}}, 
{{{fWm, FHdm}, {fWp, FHup}}, {{Lm,Um}, {Lp,Up}}}} 

The different parts are 
• Mixing of down squarks: 

jl OH* n q(I+3)i 

a L — °D u m,i, (J>Jt — >J D ^ m ,i- 



Mixing of up squarks 

^>U u m,i, U R — u ' m,ii 



~I _ qlijj . ~J _ r.( / + 3 ) i *rr* 



Mixing of sleptons 



~I _ nil* 771 ~i _ c( / + 3 )* 



• Mixing of neutral, scalar Higgs 



h J V <t>u 
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• Mixing of neutral, pseudoscalar Higgs 



Mixing of charged Higgs 



Hr\ = z ±(r d 

With cj>- = (0+)*. 

• Mixing of neutralinos: 

(B,W ,H®,H®) mix due to the matrix Z N to A° which form Majorana 
fermions 

A?" 



X 



The last entry is the mixing of the charginos: The gauge eigenstates 
W~ (fWm) and Hj (FHdm) are mixed to the negative charged mass 
eigenstates A - (Lm), while W + (fWp) and (FHup) form the new 
eigenstates A + (Lp). The new and old eigenstates are connected by 
the mixing matrices U~ (Um) and U + (Up). 

(5)-(E) 

(e) The new gauge fixing terms are 

{{Der[VP] , - 1/(2 RXi[VP])}, 
{Der [VWm] +\ [Imaginaryl] Mass[VWm] RXi [VWm] Hpm[{l}] , 

- 1/ (RXi [VWm])}, 
{Der[VZ] +Mass[VZ] RXi [VZ] Ah[{l}], 

- 1/(2 RXi[VZ])}, 
{Der[VG], - 1/(2 RXi[VG])}} 

8. No particles should be integrated out or deleted 

IntegrateOut={} ; 
DeleteParticles={} ; 

9. The Dirac spinors are the following 

dirac[[l]] = {Fd, FdL, FdR}; 
dirac [[2]] = {Fe, FeL, FeR}; 
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dirac[[3]] = {Fu, FuL, FuR}; 

dirac[[4]] = {Fv, FvL, 0}; 

dirac[[5]] = {Chi, L0, conj[L0]}; 

dirac[[6]] = {Cha, Lm, conj[Lp]}; 



A. 1.2 Parameter and Particle - File 

Additional properties and information of the parameters and particles of a model 
are saved in the files parameters .m and particles .m. A entry in the parameter file 
looks like 

{Yu, { LaTeX -> "Y~u" , 
Real -> True, 
Form -> Diagonal, 
Dependence -> None, 
Value -> None, 

LesHouches -> {{{1,1}, {Yu,l,l}}, 
{{2,2}, {Yu,2,2}}, 
{{3,3}, {Yu,3,3}}} 

}} 

and contains information about the numerical value, the position in a Les Houches 
accord file or the dependence on other parameters. Also simplifying assumptions like 
only real and diagonal values can be given to a parameter and the k)Tj^X name is 
defined. 

The particles file contains entries like 

{Su , { RParity -> -1, 

PDG -> {1000002,2000002,1000004,2000004,1000006,2000006}, 

Width -> {0, 0, 0, 0, 0, 0}, 

Mass -> Automatic, 

FeynArtsNr -> 13, 

LaTeX -> "\\tilde{u} M , 

OutputName -> "um" }}, 

and defines properties of all particles like the R-parity or the mass. 

A.2 From MSSM to NMSSM 
A. 2.1 Changing the Model File 

We want to show at one example how easy it is, to generate model files for new SUSY 
models with SARAH. For this reason we want to change the existing model file of 
the MSSM to a model file for the NMSSM (see [0| and references therein). This is 
done by: 
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1. A gauge singlet is added to the model. This is performed in SARAH by ex- 
tending the particle content with 

Sing [[4]] = {Sing, 1, S, 0, 0, 0}; 

This statements means, that a superfield called S is added to the particle 
content, which has one generation and does not couple to the three gauge 
groups. 

2. There is a superpotential term, which describes the interaction between the 
two Higgs superfields H u and H d and the gauge singlet S. Also, there is a 
cubic self coupling of the singlet field ^kS 3 . This terms are added in SARAH 
by 

TriW = {..., {Hu, Hd, S, \ [\Lambda] , 1},{S,S,S,\ [\Kappa] , 1/3}}; 

3. First, we consider a CP conserving NMSSM. This is done by declaring k and 
A as real in the parameter file: 

{\ [Kappa] , { . . . , Real -> True , . . . }} ; 

{\ [Lambda] , { . . . , Real -> True , . . . }} ; 

4. The scalar component of the singlet superfield receives a VEV after EWSB and 
splits into scalar and pseudo scalar part: 

S = -j= (05 + i<rs + vs) 
In the model file this reads as 

{SSing , {{phiS , 1/Sqrt [2] } , {sigmaS , I/Sqrt [2] } , {vS , 1/Sqrt [2] }} ; 

5. It is known, that in the CP conserving case, the scalar and pseudo scalar sector 
does decouple, i.e. the pseudo scalar singlet mixes with the pseudo scalar Higgs 
and the scalar component with the scalar Higgs: 

{{phid , phiu , phiS} , {ZH , hh}} , 
{{sigmad, sigmau, sigmaS} , {ZA, Ah}} 

The new eigenstates are CP-even h (hh) and CP-Odd A h (Ah) Higgs and the 
corresponding mixing matrices are called Z H (ZH) and Z A (ZA). 

6. The fermionic component of the singlet fields mixes with the other neutral 
fermions to five neutralinos 



10 



{{fB, fWO, FHdO, FHuO, conj [FSing] } , {ZN, Chi}} 



conj assigns the complex conjugation. We have to use it here, because we 
defined the singlet as right handed anti field. 

7. Last step is to define a Dirac spinor from the Weyl spinor for the gauge eigen- 
states 

dirac [[12]] = {FS, FSing, conj [FSing] } ; 

Calculating Vertices After evaluating the model file due to the Start-command, 
we can e.g. calculate the coupling between the scalar Higgs and the d-quarks by 

Vertex [{bar [Fd] , Fd, hh}] 

and get 

{{bar[Fd[{gtl, ctl}]], Fd[{gt2, ct2}] , hh[{gt3}]}, 
{((-I)*conj [ZH[gt3, 2]]*Delta[ctl, ct2] *Delta[gtl , gt2] * 

Yd[gt2, gtl])/Sqrt[2] , PL}, 
{((-I)*conj [ZH[gt3, 2]]*Delta[ctl, ct2] *Delta[gtl , gt2] * 

Yd[gtl, gt2])/Sqrt[2] , PR}} 

SARAH writes first a list with the involved particles and inserts the indices, which 
are used in the result: gtX for generation and ctX for color indices. The result 
consists of two entries for left (PL) and right (PR) polarization. 



Tadpole Equations SARAH also automatically calculates the minimum condi- 
tions for the vacuum: 

dV 

dvi 

We can look at the first tadpole equation after EWSB = via 
TadpoleEquationsEWSB [ [1] ] 

The output, left hand side of the equation, reads as 

mHd2*vd + (gl~2*vd~3)/8 + (g2~2*vd~3)/8 - (gl~2*vd*vu~2)/8 - 
(g2~2*vd*vu~2)/8 - (3*vS~2*vu*\ [Kappa] *\ [Lambda] )/2 + 
(vd*vS~2*\ [Lambda] "2) /2 + (vd*vu~2*\ [Lambda] ~2) /2 - 
(vS*vu*A [\ [Lambda] ] ) /Sqrt [2] 
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Ghost Interactions The neutral Goldstone boson is the first generation of the 
pseudo scalar Higgs, i.e. A\, and r] + is the corresponding ghost to W + = (W~)*. 
Thus, if we have defined the gauge fixing terms, we can also calculate ghost vertices 
like v + 7] + Al by 

Vertex [-[bar [gWmC] , gWmC , Ah [{1}] }] 
with the following result 
{{bar[gWmC] , gWmC, Ah[{l}]}, 

{-(g2*Mass[VWm]*(ZA[l, 2] *Cos [\ [Beta] ] - ZA[1, 1] *Sin [\ [Beta] ] ) )/2, 1}} 
A. 2. 2 Parametrization of the Pseudo Scalar Sector 

The first implementation of the NMSSM in SARAH runs fine, but it is also common, 
do define two rotations in the pseudo scalar sector: First isolate the Goldstone by 

a d \ / cos/3 sin/5 0\ / G \ 
a u = I — sin (3 cos (3 I I a I 
as) \ 1/ \a s ) 

and rotate afterwards (cr, erg) . 

Now we want to show how to do this in SARAH: 

1. Define a first, temporary rotation 

Rotation [[1]]={TEMP, ... 

{{sigmau, sigmad}, {AhT, ZT}}, ...} 

to temporary eigenstates A\ with a mixing matrix Z T . 

2. Parametrize this mixing matrix by 

{ZT, { Real->True, 

Dependence -> {{Cos [\ [Beta] ] , Sin [\ [Beta] ]} , 
{-Sin [\ [Beta] ] , Cos [\ [Beta] ] }} 

>} 

3. Make a second rotation to the electroweak eigenstates 

Rotation [ [2] ]={EWSB, {{AhT, sigmaS} , {Ah, ZA}>, ... } 

4. Parametrize the rotation of the second pseudo scalar Higgs and the pseudo 
scalar singlet by an mixing angle (f>: 

{ZA, Real -> True, 

Dependence -> {{1, 0, 0}, 

{0, Cos [\ [Phi]], -Sin [\ [Phi]]}, 

{0, Sin [\ [Phi]], Cos [\ [Phi]]}} }} 
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A.2.3 Adding CP Violation 

So, now we want to drop the assumption about CP conservation in the Higgs sector. 
This is done just by setting 

i\ [Kappa] , {..., Real -> False,... }}; 

{\ [Lambda] , {..., Real -> False,... }}; 

This leads to the fact that the Higgs matrix is no longer reducible and we have to 
define it as 

{{phid , phiu , phiS , sigmad , s igmau , sigmaS} , {ZH , h}} 

Analyzing Mass Matrices We want to have a look at the 6x6 mass matrix of the 
Higgs sector in SARAH. This mixing was the 4. mixing we defined in the model file, so 
we can access the corresponding mass matrix just by MassMatricesFullEWSB [ [4] ] . 
The mixing of the scalar down Higgs with the pseudo scalar up Higgs is the (1,5) 
component of this matrix. Therefore, 

MassMatricesFullEWSB [[4,1,5]] 

returns 

( (-1/2) *vS*A [\ [Lambda] ] ) /Sqrt [2] 
+ ((I/2)*vS*conj [A[\[Lambda]]])/Sqrt[2] 
- ( (3*1) /4) *vS~2*\ [Lambda] *conj [\ [Kappa] ] 
+ ( (3*1) /4)*vS"2*\ [Kappa] *conj [\ [Lambda]] 

This is the expected result: We see that for real values of k and A the contributions 
would cancel and the matrix decouples again as it should be. 

A.2.4 CKM Basis 

Although we have added CP violation in the Higgs sector, we have neglected so 
far the fact of flavor violation in the standard quark sector. But it is known that 
the Yukawa matrices are not diagonal in the gauge eigenstates of the quarks. This 
will lead after EWSB to a rotation of the gauge eigenstates to a new basis, which 
diagonalizes the Yukawa matrices: 

d L = Vj° L , u L = Vj° R 
da = U d ii° L , u R = U u u° R 

where we have assigned the gauge eigenstates with a upper index 0. We can define 
diagonal Yukawa Y u , Y d matrices in this new basis by 

Y d = U a Y°V d , Y d = UtY d °V d 
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Furthermore, in SUSY theories the soft breaking parameters are rotated in a similar 
way: 

m 2 q = Vjm° g 2 V d , m 2 u = Ulmfu u) m\ = U\mfu d 
A d = U\Afv d , A u = UlAfv u 

After this redefinition the mixing matrices V d ,V u ,U d ,U u disappear from the La- 
grangian and only the product of V d and V u stays: 

VcKM = VlV d 

Now, we want to implement this steps in our model file in SARAH: 

1. For clearance, we name the gauge eigenstates as dLO, uLO, dRO and uRO. Also 
the unrotated Yukawa matrices get an index 0. The superpotential can now be 
written as 

{{q0, Hu, uO, Yu0,l}, {q0, Hd, dO, YdO, -1},...}; 

2. Afterwards, we rotate the quarks and squarks to the CKM basis 

Rotation [ [1] ] ={SCKM, ... 

{{SdLO}, {SdL, Vd}}, {{SuLO}, {SuL, Vu}}, 

{{SdRO}, {SdR, Ud}}, {{SuRO}, {SuR, Uu}}, 

{{•[FdLO}, {FdRO}}, {{FdL,Vd}, {FdR,Ud}}}, 
{{{FuLO}, {FuRO}}, {{FuL,Vu}, {FuR,Uu}}} 

3. Now we define the rotated parameters 

{YuO, {Dependence -> 

sum [iOOl ,1,3] *sum [i002 ,1,3] *Yu [iOOl , i002] *Delta [iOOl , i002] * 
Vu [iOOl , indexl] *conj [Uu [i002 , index2] ] }} , 

{mq02, {Dependence -> 

sum [iOOl ,1,3] *sum [i002 ,1,3] *mq2 [iOOl , i002] *Delta [iOOl , i002] * 
Vd [iOOl , indexl] *conj [Vd [i002 , index2] ] }} , 

4. And the CKM matrix as product of the left mixing matrices 
{CKM, {MatrixProduct -> {Vd,Vu}}} 
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A. 3 Adding bilinear R-Parity Violation 

As next step we want to consider the possibility of bilinear R- Parity violation. This 
can be done in the following way: We write a superpotential coupling between the 
gauge singlet, the left handed leptons and the up Higgs. If we now consider the gauge 
singlet as right neutrino, we get the /if SSM [[H], [Tj|. 

• Define the right handed neutrinos 
Sing [[4]] = {vR, 3, v, 0, 0, 0}; 

• Define the superpotential 

TriW = {. . . ,{1, Hu, v, Yv, 1}, 

{v ,v, v,\ [Kappa] ,1/3}, 
{Hu, Hd, v, \ [Lambda] ,-1}}; 

• Give the left and right handed neutrinos a VEV 

{SvL, {vL, 1/Sqrt[2]}, {sigmaL, I/Sqrt [2] } , 

{phiL, 1/Sqrt[2]}}, 

{SvR, {vR, 1/Sqrt[2]}, {sigmaR, I/Sqrt [2] } , 

{phiR, 1/Sqrt[2]}} 

• Define the mixings in the Higgs sector. The sneutrinos mix with the neutral 
Higgs and the selectrons with the charged Higgs: 

{{phiu, phid, phiR, phiL}, {hh, ZH}}, 
{{sigmau, sigmad, sigmaR, sigmaL} , {Ah, ZA}}, 
{{conj [SHup] ,SHdm, SeL, SeR},{Hpm,ZP}} 

• Define the mixings in the fermion sector. The neutrinos mix with the neutrali- 
nos and the eletrons with the charginos: 

{{fB, fWO, FHdO, FHuO, conj [FvR] , FvL}, {L0, ZN}}, 

{{{fWm, FHdm, FeL}, {fWp, FHup, conj[FeR]}}, {{Lm,Um}, {Lp,Up}}} 

A. 4 Adding an additional Gauge Group 

Let's say, we want now to add an additional U(l) to the model. This is done in 
SARAH just by adding a new gauge superfield: 

Gauge[[4]]={X, 1, extra, g4}; 
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The name of the superfield is just X, so the new vector boson is called VX and the new 
gaugino f X. The dimension of the group is 1 and we called it extra in comparison 
to hypercharge or color. The parameter for the coupling constant is g4. 
As second step we have to define the charge of all chiral under this new group,e.g. 

Doub[[4]] = {Hup, HuO, 1, Hu, 1/2, 1, 0, -1/2}; 

This means up Higgs has hypercharge |, transforms under SU(2) L , does not trans- 
form under SU(3)c and has charge — \ with respect to the new gauge group. 
Obviously, the Higgs couple to the new gaugino X (fX in SARAH). This leads to a 
further expansion of the neutralino mass matrix by 

{{fB, fWO, FHdO, FHuO, conj [FvR] , FvL, fX}, {L0, ZN}} 

But, when we start now the evalution in SARAH, we get this messages: 

WARNING! U(l)~3 (extra) Gauge Anomaly! 

WARNING! U(1)*SU(2)~2 (extra x left~2) Gauge Anomaly! 

That means, we have not been very careful in choosing the new gauge quantum 
numbers and should improve this choice. 
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